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ABSTRACT 



Aims. We analyze observational data from 4 instruments to study the correlations between chromospheric emission, spanning the 
heights from the temperature minimum region to the middle chromosphere, and photospheric magnetic field. 

Methods. The data consist of radio images at 3.5 mm from the Berkeley-Illinois-Maryland Array (BIMA), UV images at 1600 A 
from TRACE, Can K-line filtergrams from BBSO, and MDI/SOHO longitudinal photospheric magnetograms. For the first time 
interferometric millimeter data with the highest currently available resolution are included in such an analysis. We determine various 
parameters of the intensity maps and correlate the intensities with each other and with the magnetic field. 

Results. The chromospheric diagnostics studied here show a pronounced similarity in their brightness structures and map out the 
underlying photospheric magnetic field relatively well. We find a power law to be a good representation of the relationship between 
photospheric magnetic field and emission from chromospheric diagnostics at all wavelengths. The dependence of chromospheric 
brightness on magnetic field is found to be different for network and internetwork regions. 
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1. INTRODUCTION 

The solar chromosphere is known to be highly structured at 
different spatial scales and the range of chromospheric inho- 
mogeneities is best revealed when observed at different wave- 
lengt hs from the ex treme ultraviolet to millimeter waves (see, 
e.g., Solanki 2004). Even far from activity complexes, in the 
quiet Sun magnetic field is generally believed to be responsi- 
ble for chromospheric structuring. This is consistent with the 
strong correlation found between excess emission in the cores 
of the Can K and H resonance lines (the primary chromo- 
spheric diagnostic) and the presence a nd strength of quiet - 
Sun magnetic flux (e.g. iLeightonl [19591: ISchriiver et all fl989). 
Consequently, these resonance lines commonly serve as indica- 
tors of changes in the chromospheric structure related to global 
magnetic activity a nd heating of the outer atmosphere of the 
Sun a nd cool stars dRutten & Uitenbroeklll99lt ISchriiver et alj 
1989). In additi on to the c hromospheri c emission related to the 
magnetic field, ISchriiven 1 (119871 [[992) has suggested that part 
of the radiative flux from solar and stellar chromospheres, re- 
ferred to as basal flux, has a nonmagnetic origin and is believed 
to stem from acoustic hea ting of the chromosphere (see, e.g. 
iRutten & Uitenbroeklll99ll) . 

The primary brightness patterns revealed in two-dimensional 
images obtained in the ionized calcium H and K lines include 
bright plages, chromospheric network located at the boundaries 
of supergranular cells, and bright points in the interior of the 
cells - internetwork grains or K2V bright points. Chromospheric 
network features are known to display a one-to-one spatial cor- 
respondence with regions of enhanced photospheric magnetic 
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field (Skumanich et al. 1975; Sivaraman & Livingston! Il982b 
iNindos & Zirinl [T998I : ISivaraman e t al. 2000), while the (mag- 
netic) origin of the internetwork grains is a controversial is- 
sue. From the analysis of simultaneous Kitt Peak Can K- 
line spectroheliograms and photospheric magnetic area scans, 
ISivaraman et alj 12000) showed the association of internetwork 
K2V bright p oints with magnetic featu r es, co nfirming the ear- 
lier results of Sivaraman & Livingston (1982). However, from 
analysis of Ca n K-line and magnetic field data from the BBSO, 
INindos & Zirinl d!998l) have found that not all internetwork 
bright points are associated with magnetic fields. The existence 
of 2 distinct classes of cell-interior grains, which differ in tempo- 
ral behaviour and spatial motion and have a different origin, was 
revealed by observations wit h the Swe dish Solar Telesco pe re- 
ported by Bra ndt et al.1 dl992l) , whereas lLites et al.l (11999) found 
no obvious rela tion of the internetwo rk brightenings to the mag- 
netic field, and Remling et al.l d 19961) concluded that they are a 
totally nonmagnetic phenomenon. 

Quantitative studies of the relationship between magnetic 
field and the chromosphere (represented by the calcium K and 
H-lines), in both the active and quiet sola r atmosphere, have been 
carried out by a number of authors (e.g . Skumanich et al. 19751: 
I Schriiver etakl 119891: INindos & ZMnMl998l: jHarvev & White! 
Il999h lOrtiz & Rastl [2005t iRezaei et all 120071) In the quiet- 
Sun ch romosphere[Skumanich et alj d 19751) and INindos & Zirinl 
d!998l) found a linear correlation between the K-line inten- 
sities and the absolu te values of the magnetic field strength. 
ISchriiver et al.l d 1989b extended this relationship to solar active 
regions and suggested that the relation between the K-line ex- 
cess flux density (after the subtraction of the basal flux from 
the observed Ca n K-line flux) and the magnetic field flux den- 
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sity is best described by a power law with an exponent of about 
0.6. Considering 4 types of brightness structures, which included 
active regions, decaying a ctive regions, the enhan ced network 
and quiet (weak) network, Harvey & White] d 19991) ascertained 
that in both active and quiet Sun the relation between the mag- 
netic flux density and calciu m residual intensity f ollows a power 
law with an exponent of 0.5. lQrtiz & Rasj d2005l) c ited a power- 
law ex ponent value of 0.66 for the quiet Sun, while iRezaei et al.1 
(120071) obtained a lower exponent of 0.2 for locations in a quiet- 
Sun region, and higher values of 0.4-0.5 for network locations, 
de pending on th e magn etic field threshold value. 

RuttenetdJdl991) found power-law relations between ex- 
cess flux densities in stellar chromospheric lines with differ- 
ing temperature responses, with the exponents increasing with 
increasing difference between their temperatures of formation. 
Relations between fluxes originating in the same atmospheric 
regime were found to be close to linear (i.e., having a power- 
law exponent of unity). For the quiet Sun the high degree of co- 
spatiality of the brightness features seen in the Ca n K line and 
the TR ACE UV channels has been demonstrated bv lRutten et al.l 
(119991) . 

Almost all these studies have been limited to the vis- 
ible and UV spec t ral ran ge. An exception is the work of 
iLindsev & Jefferiesl dl99ll) . who studied the cross-correlation 
statistics of the brightness variations in the submm (850//m) and 
Can K-line images. Their results are confined to the tempera- 
ture minimum region and the heights of the low chromosphere. 
No detailed studies have compared millimeter emission, which 
originates from the middle chromosphere, with other chromo- 
spheric indicators. From observations at 350, 850 and 1200 
yum with t he James Clerk Ma xwell Telescope at a resolution of 
about 20". lLindsev et al.ldl995l) reported submillimeter chromo- 
spheric brightness structures which resemble the supergranular 
network, although the resolution was marginal for such a study. 
Initial qualitative results of the comparison of chromospheric 
structure, identified in millimeter images obtained with BIMA 
at 3.5 mm with a resolution of about 10", with calcium images 
from BBSO and line-of-si ght photospheric m agnetograms from 
SOHO/MDI were given in I White et all d2006l) . 

In this paper we present an observational study of the rela- 
tionships between the photospheric magnetic field and emission 
at different temperatures from chromospheric heights. We study 
the chromospheric observational data referring to the tempera- 
ture minimum region (UV continuum at 1600 A), low chromo- 
sphere (Ca n K line) and middle chromosphere (continuum at 3.5 
mm). This is the first analysis of this type that includes millime- 
ter observational data. Unfortunately, simultaneously recorded 
observational data on the chromospheric magnetic field suitable 
for such an analysis are unavailable. For this reason, we focus on 
the MDI/SOHO longitudinal photospheric magnetograms and 
study the dependence of the chromospheric intensity at differ- 
ent wavelengths on the photospheric magnetic field strength. 

In Sec. 2 and 3 we describe the data and instruments used 
for the analysis and provide the details of the data reduction. 
In Sect. 4 the topography of the quiet chromosphere, chromo- 
spheric intensity distributions and cross-correlations between the 
magnetic field and chromospheric emission are presented. In 
Sect. 5 we summarize and discuss the results. 

2. OBSERVATIONAL DATA 

We analyze observations from the Berkeley-Illinois-Maryland 
Array (BIMA) at 85 GHz (A = 3.5 mm), Can K-line fil- 
tergrams from the Big Bear Solar Observatory (BBSO), UV 



images at 1600 A from the Transit ion Region and Coronal 
Explorer (TRACE, lHandv et al.ll 19991) and magnetograms f rom 
the Michelson Doppler Imager (MDI, IScherrer et al.l 1 19951) on 
the Solar and Heliospheric Observatory (SOHO) satellite, ob- 
tained on May 18, 2004. The data consist of temporally resolved 
images of the measured intensities, which are given in Kelvin for 
millimeter data (see Sec. 12. It , in relative units for Can K-line 
and 1600 A data (Sec. 12.21 and 12.31 respectively), and in Gauss 
for MDI photospheric magnetograms (Sec. 12.41 ). The overview 
of the data (including instruments, operating wavelengths, spa- 
tial resolution, cadence and image pixel size) used in the present 
study is given in TableQ] 



2. 1. BIMA observations at 3.5 mm 

The 3.5 mm observations of a quiet-Sun area near the center of 
the disk were acquired with 9 BIMA antennae during the time 
interval 16:43-20:21 UT (~3.5 hours) when atmospheric phase 
conditions were good. The data were taken at a cadence of 10 
seconds. The interferometer was in its most compact ("D") con- 
figuration for these observations, providing a spatial resolution 
of order 12" within a 2' field of view (FOV). The images were 
deconvolved using a maximum entropy method (MEM) and re- 
stored with a Gaussian beam. The details of the interferomet- 
ric image restoration of the BIMA data, including description of 
th e Umitations impo sed by the MEM deconvolution, are given 
in I White et al.l d2006l) . The final images were 128 pixels square, 
with a cell size of 3" and the flux in the inner 80 pixels of the 
FOV with the center at (-17", 71") in heliocentric coordinates. 
All images were corrected for the primary beam response and 
flux units were converted to brightness temperatures for a 12" 
beam. Data within a 96" radius of the pointing center are used 
in the analysis below. The zero level in the BIMA image corre- 
sponds the coldest feature on the surface in the chromosphere 
that is optically thick at 3.5 mm wavelength. This should be at a 
temperature of order 6500-7000 K. The temperature fluctuations 
in the BIMA image represent direct temperature measurements 
above this minimum value on the optically thick surface in the 
chromosphere, ranging from to 1 100 K in the map for the en- 
tire 3.5 hour period. 



2.2. Ca n K line observations at 393.3 nm 

The Can K-line observations were obtained with the 25-cm re- 
fractor of BBSO. A 5h sequence of filtergrams in the core of the 
Ca II K-line (393.3 nm) with a bandpass of 0.6 A was recorded 
by a KX-260 CCD camera, which has 512 x 512 pixels and a 
dynamic range of 14 bits. The exposure time was 250 msec, 
the cadence was 30 sec, recorded field of view was 300" x 
300" at a pixel size of 0.645". The observations were taken 
under moderate to fair seeing conditions. The K-line data were 
preproces sed, which includ ed calibration by dark and flat-field 
frames (iDenker et al.l l200lT) . The intensity scale of the images 
was equalized by setting the quiet-Sun integrated background 
intensity in all images to the same value. No absolute flux cali- 
bration was available for the data, hereafter we use the Ca n K- 
line intensities in arbitrary units. We note that the bandwidth of 
the calcium filter tuned to the line center was too large to sep- 
arate the K-line peaks and thus the images obtained depict the 
structure of the solar chromosphere collected from a wide range 
of heights, but dominated by mid-chromospheric contributions, 
as sampled by the K2 peaks. 
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Instrument 


Passband 


A 


AA(A) 


Spatial resolution (") 


Time cadence 


Pixel size (") 


MDI/SOHO 


Nil 


6768A 


0.094 


4 


1 min 


1.96 


TRACE 


C IV+cont 


1600 A 


275 


1 


30 min 


0.5 


BBSO 


CanK 


3933A 


0.6 


2-4 


30 s 


0.645 


BIMA 


continuum 


3.5 mm 




12 


10 s 


3.0 




Fig. 1. Portrait of the solar chromosphere averaged over 3.5 h at the center of the Sun's disk at 4 different wavelengths on May 
18, 2004. From top left to bottom right: MDI longitudinal photospheric magnetogram (resolution of 4"), TRACE 1600 A image 
(resolution of 1"), Can K-line center image from BBSO (resolution of 2-4") and BIMA image at 3.5 mm (resolution of 12"). The 
images are created by averaging over the 3.5h period of the BIMA observation. The details of the data units are given in Sect. [2] 
Dashed circles mark the 96" radius BIMA FOV where we regard fluxes to be reliable. X and Y axes are in arcsec from the disk 
center, which is located at (0, 0). 



2.3. TRACE observations at 1600 A 

In parallel with the BIMA observations a total of 14 TRACE 
images with a cadence of 20 minutes were recorded at a wave- 
length of 1600 A at full resolution (0.5" pixel size). The im- 
ages were preprocessed, which included subtraction of the dark 
pedestal and current, as well as division by a corrected flat 
field. The data were normalized to the exposure time, giving 
data numbers per second (DNs~ l ) with no absolute calibra- 
tion. Both standard models of the quiet-Sun atmosphere such 



as FAL C dFontenla et alj 1 1993b and no n-LTE radiation hy- 
drodynamic simulations such as those of Fossum & Carlsson 
d2005l) place the disk-center radiation at 1600 A in the tem- 
perature minimum reg ion at a height of around 500 km 
above the photosphere dHandvetaT]|l999l:lKriigeretalJl200ll: 
iFossum & Carlssonll2005b . although a contribution from 10 5 K 
transition-region gas is also present due to the C IV lines. In the 
quiet Sun this contributio n is expected to be small (lHandv et alj 
ll998HKriigeretal.ll2001l) . 
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Fig. 2. The same as Fig.[T]plotted at the BIMA resolution of 12". In the top left panel the absolute value of the MDI magnetogram 
is depicted. Overlays of contours of 3.5 mm emission at (30, 50, 70, 90, 95, 98)% of maximum brightness are overplotted in white. 



2.4. MDI magnetograms 

One-minute ca dence full-disk MDI magnetogramsQ (spatial res- 
olution of 4", IScherrer et alHl995h were selected for the 3.5 h 
BIMA observing period, which resulted in 188 images of the 
proxies of the photospheric magnetic field. The noise in the 
original 1-min MDI magnetograms is typically about 20 Gauss 
(Sche rrer etafll 19951) . which is comparable to the magnetic sig- 
nal of quiet-Sun network features. By averaging in time over 
24 consecutive individual magnetograms (30-min sequence) we 
reduced the noise to 4 Gauss. Averaging over the whole 3.5 h 
observing period (188 images) resulted in a level of the random 
noise of about 1.5 Gauss. Prior to averaging the influence of so- 
lar rotation on the time series was removed (see Sect. [3). 



3. DATA COALIGNMENT AND REDUCTION 

We formed 4 observational data cubes I(x, y, ?), which were sub- 
sequently aligned using a two-step procedure. The first step in- 
volved correction of differential solar rotation by using standard 
procedures of SolarSoftware in the IDL environment. All im- 
ages were rotated to the fiducial time of 20:00 UT. In the second 



1 lev 1.8 data. 



step the shift of the images was determined and applied by a 2D 
cross-correlation of each image of the sequence to the first image 
of the sequence. 

After a sub-arcsecond accuracy was reached in the align- 
ment of the data cubes we averaged the images in time over the 
whole observing period (16:43-20:20 UT) to obtain the mean 
(time-averaged) MDI magnetogram, Ca n K filtergram, TRACE 
1600 A image and BIMA image at 3.5 mm. Additionally we av- 
eraged over the uninterrupted sequence of images of good qual- 
ity within the time interval 19:00-19:30 UT to check the stability 
of the large-scale quiet-Sun structures. We based our correlation 
analysis on the mean images to avoid the scatter introduced by 
the chromospheric oscillations within cell interiors and the chro- 
mospheric network. Consequently, by averaging the images we 
also suppress proper motions of the K-grains and of bright points 
in the network. 

Finally, the full disk B\\ MDI images were chosen as the co- 
ordinate reference against which all other mean images were co- 
aligned by using cross-correlation methods. Absolute values of 
the MDI magnetograms were formed prior to alignment. The re- 
maining image displacements do not exceed 1" and we believe 
that they do not affect the results of our analysis. A 3.5h time- 
averaged image of the BIMA FOV made at each of the 4 wave- 
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length bands is shown in Fig.Q]at the original spatial resolution 
of each band. 

The cross-correlation analysis between the chromospheric 
images and photospheric magnetograms was done both for the 
MDI magnetogram resolution (4") and for the resolution of the 
BIMA data(12"). In the former case the Ca n K-line and TRACE 
1600 A images were rebinned to make the size of the pixels 
equal to that of the magnetograms. In the latter case we degraded 
all images to the resolution of the BIMA data prior to cross- 
correlating them. The shift between the BIMA images and the 
MDI magnetograms did not exceed 1 BIMA image pixel (< 3"). 
Figure [2] shows example time-averaged images at BIMA resolu- 
tion. 



4. RESULTS 

4.1. Mapping of chromospheric structure 

FigureQ]depicts the time-averaged quiet solar chromosphere at 3 
different heights from the temperature minimum region (TRACE 
1600 A image) through the lower chromosphere (Can K-line 
filtergram) to the middle chromosphere (BIMA image at 3.5 
mm) at the original resolution of the instruments (see Table [TJ 
as well as the corresponding MDI/SOHO line-of-sight photo- 
spheric magnetogram. Since the QS region being analyzed lies 
near the disk center, the line-of-sight magnetic field, measured 
by MDI, can accurately be regarded as the vertical component 
of the magnetic field vector. Of particular interest is the relation 
of the BIMA image to the others. Visual inspection of Figure- 
veals strong morphological similarities between different chro- 
mospheric emissions and the photospheric magnetic field, most 
prominently on spatial scales of the enhanced magnetic flux and 
chromospheric network (brightenings in all 4 frames). In snap- 
shot images (not plotted) brightenings appear not only at su- 
pergranular cell boundaries, but also within the cells. Averaged 
in time, bright network becomes more pronounced, while low- 
intensity short-lived internetwork bright points are significantly 
suppressed and eventually vanish with increased averaging. In 
general, use of time averages reduces differences between the 4 
images due to the brightness variability of the structures. 

The impression of a high degree of correlation of photo- 
spheric magnetograms (absolute value) with the Can K-line, 
1600 A and 3.5 mm images is confirmed quantitatively. Within 
a radius of 96" the linear correlation coefficient between the ab- 
solute magnetogram and calcium image is 0.75, while it is 0.63 
between the magnetogram and emission at 1600 A. 

The degree of spatial agreement between the radio images 
and the other chromospheric diagnostics and the photospheric 
magnetic field is more clearly visible in maps degraded to the 
BIMA resolution of 12", shown in Fig. [2] with radio contours 
overlaid on the other images at (30, 50, 70, 90, 95, 98)% of 
maximum mm brightness. The correlation coefficients between 
the mm image and the magnetogram are 0.56 (linear) and 0.72 
(rank), while for the mm image and Ca n K-line image they are 
0.65 and 0.69, respectively. At a resolution of 12" the magne- 
togram shows a higher correlation with the calcium and 1600 A 
images than at 4" resolution, with the linear correlation coeffi- 
cient being 0.82 and 0.74, respectively. 

4.2. Histograms of chromospheric intensity distribution 

With sufficiently good statistics the brightness histogram is ex- 
pected to have a log-normal shape, comparable with a single 
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Fig. 3. Intensity distribution histograms over the 2-D quiet-Sun 
region near the disk center obtained, from top to bottom, in the 
TRACE 1600 A channel, the Can K line and at 3.5 mm. In all 
three panels dotted curves depict log-normal fits. 



population dPauluhn et al.l H000). To test whether this is the 
case for the chromospheric intensities studied here we built in- 
tensity distribution histograms (observed frequency distribution 
normalized to unit area) for the selected 2D images at 1600 A 
and 3.5 mm and fitted t hem with a single l og-normal distribu- 
tion function. Following Fontenla et al.l d2007l) histograms of the 
Ca n K intensity were fitted with a modified log-normal function 
having an additional parameter - an offset which represents line 
scattering of the photospheric background radiation. In Fig.|3]we 
show the derived histograms for the observed intensities at 1600 
A, in the Ca n K line and at 3.5 mm and their log-normal fits with 
solid and dotted curves respectively. 

A log-normal with an offset provides a good fit to the ob- 
served Ca n K line intensity distribution, although the observed 
tail enhancement is slightly higher than that of the fitted log- 
normal (see central panel in Figj3}- F° r the 1600 A intensities 
the quality of the fit diminishes at low and high flux values and 
the departure from a log-normal form increases. Finally, due to a 
substantially higher noise level in the millimeter brightness data 
(and the smaller number of resolution elements available for the 
histogram) it is difficult to evaluate the quality of the log-normal 
fit for this dataset. On the whole, a tendency for bigger depar- 
tures from a log-normal can be observed when spatial resolution 
is decreased. However, we cannot rule out the possibility that 
we are seeing the effects of different heating mechanisms in cell 
interiors and network lanes being manifested in the intensity dis- 
tributions. 



6 



Loukitcheva et al.: Chromospheric structure from multi-band observations 




Fig. 4. The Can K-line intensity vs. the absolute value of the magnetic flux density for a 2" pixel size. The images are averaged in 
time over 3.5h. The vertical line at 1 .5 G in the central panel represents the lcr noise limit of the magnetogram signal. The fit results 
are represented by a solid curve for a power-law fit with free exponent, dotted curve for a power-law fit with the exponent fixed at 0.6 
and dashed for a linear fit. The light grey solid curve outlines the Ca n K intensity binned in 70 intervals containing an equal number 
of points. Histograms on the left show the distribution of magnetic field strength in 4 subranges of the Can K-line flux. Histograms 
on top of the central panel depict the distribution of calcium intensity in 4 subranges of the magnetic flux. For each histogram the 
corresponding range of values considered is marked in the main figure: by horizontal solid lines for magnetic flux distributions and 
by vertical dotted lines for the Can K-line flux histograms. The numbers of counts, n, represented in the histograms are indicated. 
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Fig. 5. The same as Fig.g]for the TRACE 1600 A flux at the 4" resolution of MDI. 
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4.3. Relation between chromospheric emission and 
photospheric magnetic field 

In order to study the relationship between the intensity of chro- 
mospheric bright features and the intensity of the associated 
magnetic elements, we computed scatter plots for the pairs of 
images. Firstly, we considered the pointwise (pixel-to-pixel) cor- 
relation globally in the field of view without regard to the lo- 
cation of the points. In Fig. |4] we show the resulting cross- 
correlation statistics of the brightness features in the time- 
averaged calcium image and the photospheric longitudinal mag- 
netogram. The pixel intensity correlation is plotted in the main 
frame of the figure, with histograms of the calcium intensity in 4 
subranges of the magnetic flux density superposed in the upper 
frame and of the magnetic field in 4 subranges of the calcium in- 
tensity to the right. We employed averaging over time to reduce 
the noise level in the scatter plots. 

We searched for a quantitative relation between the quiet- 
Sun mean magnetic flux density (its absolute value) per pixel 
of 2" and Can K-line intensity, looking for the best least- 
squares fit among functions of 2 type s: linear (i.e. a power 
law with unit exponen t, Eq.Q] following Skuma nTch et alj|1975t 
iNindos & Zirinl 19981) and a "generalized" power law (Eq.[2] fol- 
lowing lSchriiver et al.|[l989l) . 

I = I +a-\B\ (1) 

/ = I + a ■ \B\ b (2) 

The item 7p in Eqs. [T] and [2] denotes "zero" or "basal" flux, 
in troduced by |Schrii ve J (1 1 987lh For comparison with the results 
of lSchriiveret all (11989b and of others, who obtained a power 
law exponent of about 0.6 for active and quiet regions, we also 
calculated the best fit with the power-law exponent, b, fixed at 
0.6. Fits to the data points and to binned values (averaged over 



all pixels within a narrow range of magnetic flux values) points 
(see central panel in Fig. |4j gave very similar results. Only data 
points with an absolute value of the magnetic flux exceeding the 
lcr-noise level were considered for fitting. 

In the same manner we studied the relationship between the 
observed UV intensity at 1600 A and absolute magnetic flux 
density per 2" aperture, which is depicted in Fig. [5] and between 
intensity at 3.5 mm and absolute magnetic flux density per 12" 
aperture, plotted in Fig. [6] The decrease in range of < B > is due 
to the lower spatial resolution of mm data. In Table [2] we list the 
results of the function fitting to the data points for all 3 relation- 
ships investigated, and in particular we compare y 2 values for 3 
different fitting functions obtained for Poisson weighting and list 
the power law exponents obtained. 

As seen from the scatter plots in Figs. [4] [5] and [6] all three 
quantities, Can K-line, UV at 1600 A and 3.5 mm emission, 
show a prominent non-linear dependence on photospheric mag- 
netogram signal. This conclusion is confirmed by the best-fit 
curves plotted in Figs. [4] and [6] and by the x 2 values listed 
in Table [2] which reach the smallest values for the fit in the 
form of a power law. The best-fit value of the power-law ex- 
ponent is 0.31 for the relationship between the calcium inten- 
sity and the magnetic flux density, whic h lies between t h e high 
val ue of 0.6 + . 1 and 0.66 found by ISchriiver et all ([l989) 
and lOrtiz' & Rast (2005) and the lower value of 0.2 obtained by 
iRezaei et alj d2007l) . For 1600 A intensity vs. magnetic flux we 
obtain an exponent of 0.36. These results do not change if we 
use, e.g., 10-minute or 30-minute averages instead of the 3.5h 
average used above. However, the scatter of the points at a given 
< B > decreases significantly with increasing averaging. 

The relationships and density distributions that we determine 
nicely demonstrate that most of the data points in the quiet-Sun 
region are characterized by low magnetic flux and low chromo- 
spheric emission. In other words, the majority of data points lies 
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Table 2. Least-squares deviation,^- 2 , of the best fits to the relationships between the photospheric magnetic field and chromospheric 
emissions 



Flux-flux 


Linear" 


Power law b with b = 0.6 


Power law b 






X 


x 2 


2 

X 




b 


\B\ - IcallK 


11.97 


9.17 


8.44 


0.31 ± 


0.01 


\B\ ~ AeooA 


11.53 


10.64 


10.47 


0.36 ± 


0.01 


\B\ - /3.5mm 


161.68 


132.78 


111.06 


0.20 c ± 


0.01 



a seeEq.0 
b see Eq.|2] 

c calculated for \B\ > 4 G 



in the lower left quarter of the scatter plot and "anchors" the re- 
lation. There are too few data points with high chromospheric 
emission and enhanced magnetic flux (greater than 100 G) to 
very reliably determine the form of the relation (see Figs.|4]and 
0. The rather small lcr error bars of the fitted exponents (last 
column of Table 2) are the formal errors returned by the fitting 
procedure and appear rather optimistic to us. 

Note that fitting of the data directly with a power law con- 
taining a zero flux item 7q, corresponding to emission from loca- 
tions of zero magnetic flux (non-magnetic component) is simi- 
lar to determining the zero po int of chromospheric e mission us- 
ing the algorithm applied by Schriive r et al 1 (119891) . according 
to which a range of test values of Iq is subtracted from the ob- 
served flux / until the mean relationship Al = / - Iq - a ■ \B\ h 
reaches its highest correlation by minimizing the perpendicular 
offsets between the excess flux and the fit. We have tried fitted 
the data using both algorithms and obtained similar zero flux val- 
ues and power law expon ents. To demons t rate th e success of the 
algorithm, introduced by Schriiv er et al.l {T989), applied to the 
data analyzed here, we plot in Fig.|7]the excess Can K-line and 
1600 A fluxes after subtraction of the derived basal flux (listed 
in the Y-axis captions) and the corresponding power-law curves. 
These fits yield zero-flux values Iq and power law exponents b 
very similar to those derived when fitting the observed fluxes 
directly with a function of the form I = Iq + a ■ \B\ b (listed in 
Table©. 

The plot of millimeter brightness vs. magnetic field strength 
displays a considerable scatter about the mean relation, which 
may have different causes (see Fig. [6]). Pixels in the lower left 
corner with very low mm brightness can result from the under- 
estimation of the mm flux due to the high noise in the snapshot 
images (White et al. 2006). Some of the peculiar pixels with in- 
tense mm brightness (in the upper left corner) are located above 
localized (small-scale) photospheric magnetic elements and near 
polarity changes. Therefore the corresponding magnetic flux al- 
most vanishes when degrading the magnetogram to the BIMA 
resolution. Again, a power law gives a far better representation 
of the data than the linear curve. In this case, a rather low expo- 
nent of 0.20 is obtained. 

4.4. Magnetogram signal thresholding experiment 

In Figs.|U|5]and|6]the brightness shows signs of increasing more 
rapidly for smaller fields with B < 30 G. This may be saying that 
the chromospheric emission from weak-field regions (internet- 
work) has a different source than from stronger network fields. 
To study how the correlation depends on the magnetogram signal 
we have performed a thresholding experiment. We have succes- 




° -4 -2 2 4 6 
log(|B|,G) 



Fig. 7. Log-log scatter plots for a) the Ca n K-line excess mean 
intensity (after the subtraction of the zero flux) and b) TRACE 
1600 A excess mean intensity vs. the absolute value of the mag- 
netic flux density (for 2" pixel size). The solid lines represent 
the least-squares power-law fits to the data. 



sively fitted the relations between the chromospheric radiation 
and photospheric magnetic flux exceeding a number of thresh- 
old values (from to 35 Gauss for correlation at the resolution 
of the photospheric magnetogram and from to 20 Gauss for 
the resolution of the mm images). In Figs. [8]and[9]we compare 
the resuting x 1 values for linear and power law fits as a func- 
tion of the B threshold value, plotted with triangle and diamond 
symbols, respectively. We also show the dependence of the fitted 
power-law exponent on the threshold value. 

The power-law exponent of the fit to the calcium intensity 
grows with the imposed B threshold (see Fig. |S). The value of 
the power-law exponent of ^ 0.5 is reached both when we dis- 
card pixels with weak fields and consider only those pixels with 
magnetic flux exceeding 20 G, and when we include all pixels 
in the fit (including those dominated by noise). The dependence 
of the power-law index on B threshold for the intensity at 1600 A 
shows a different behaviour. It increases until it reaches a value 
of about 0.6. For magnetic flux density exceeding 13 G it starts 
to decrease rapidly for higher thresholds (Fig. |9j3). The scatter 
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Fig. 8. Thresholding experiment for the Call K-line flux vs. ab- 
solute value of the magnetic field: a) x 1 an d b) power law expo- 
nent as a function of the imposed B threshold. The vertical line 
represents the lcr-noise level of the magnetogram signal. 




l=IO+a*B A b 
l=IO+a*B"0.6 
l=a*B + b 



5 10 15 
B thres, G 



20 




5 10 15 
B thres, G 



Fig. 9. The same as in Fig.[7]for the TRACE 1600 A flux. 



of the points is so large that for these data it is not possible to 
clearly distinguish between linear and power-law fits. The sen- 
sitivity of the power law parameters to the B threshold is due to 
the small number of points at large < B > (see Fig. f5j, which is 
due to the fact that the observations refer to the quiet Sun. 

The scatter in the intensity at 3.5 mm is so large that its re- 
lationship with the magnetic flux density can be equally well fit 
by a linear or a power-law function. Also, the x 1 values of the 
fits are large. However, according to the thresholding experiment 
for the correlation between the mm intensity and the magnetic 
flux (not plotted), a power law adequately fits the data for mag- 
netic fields stronger than 3 G (~ lcr-noise level) and weaker than 
15 G. For weak fields, the large scatter affects significantly the 
quality of the fit, while for fields stronger than 15 G there are 
insufficient data pixels for confident power-law fitting. 



4.5. Discrimination between IN and NW emission. Separate 
flux analysis 

Next we study the relations between the chromospheric emis- 
sion and the magnetic flux separately for the areas represent- 
ing cell interiors or internetwork (IN) and chromospheric net- 
work (NW). In other words, we tried to separate magnetically 
weak-field (IN locations) and strong-field (NW locations) com- 
ponents of the emission. To discriminate between IN and NW lo- 
cations we employ ed a brightness threshold technique (see, e.g. 
Krijg er et al.ll200l|) . To this end we fitted the Can K-line inten- 
sity histogram constructed for the time-averaged intensity map 
by a double-Gaussian function, assuming that the intensity his- 
togram is dominated by two main populations representing the 
IN and NW intensities (see Fig.lTOl. The resulting maxima of the 
single-Gaussian components were used as intensity thresholds to 
define respectively internetwork and network areas on the mean 
map. The resulting IN (black) and NW (white) binary masks are 
shown in the bottom panel of Fig. |T0j The points belonging to 
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Fig. 10. IN/NW definition. Upper panel: Call K-line mean in- 
tensity image. Central panel: mean intensity histogram (thick 
solid), single-Gaussian profiles (dotted, dashed) from a double- 
Gaussian fit to the data points with vertical lines of the same type 
marking the location of the Gaussian maxima. Bottom panel: 
Masks used for the definition of the internetwork (black) and 
network (white) areas. 



neither of these categories are indicated by an intermediate shad- 
ing. 

We apply these binary masks to study the relationship be- 
tween the photospheric magnetic field and IN/NW fluxes, re- 
spectively (see Figs. II 11121 and [TBI . Following the approach 
outlined in Sect. 4.3 we fitted the scatter plots with linear and 
power-law functions. For the areas selected to represent the net- 
work a power law provides the best fit to the data. In Table [3] 
we list the exponents from the best-fit power laws for the im- 
ages at the resolution of the MDI magnetograms (4") and of the 
BIMA images (12"). Below we discuss the results obtained for 
the 3.5h image averages, but for comparison we list in Table [3] 
also those for the 0.5h time-averages. For the calcium NW in- 
tensity, averaged in time over 3.5 hours, the best-fit power-law 
exponent is 0.6, whic h is consistent with the value derived by 
ISchriiver et ail d!989l) for the Ca n K-line flux in an active region 
and by Harvey & White! (1 19991) for the enhanced network. We 
note a big difference from the value we derived with no image 
segmentation (0.6 vs. 0.31, see Tabled, which means a stronger 
dependence of the network calcium intensity on the magnetic 
flux. At a resolution of 12" the best fit for the NW Can K-line 
flux is also provided by a power law. This time the power law 
exponent is 0-54, whi ch is also in agreement with the results of 
Schriiver et al.l d!989l) . who state that the mean dependence of 
the Ca ii K-line excess flux on the magnetic flux density is more 
or less scale independent. Analysis of the network emission at 
1600 A gives similar results. The power-law exponent in the re- 
lationship between the NW intensity at 1600 A and the magnetic 
flux density is =* 0.5 at 4" resolution and 0.66 at 12" resolution 
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Fig. 11. Comparison of IN and NW emission averaged over 0.5h: 
histograms of a) IN and c) NW magnetogram signal, b, d) Ca n 
K-line flux vs. absolute value of the IN and NW magnetic field, 
respectively. The solid vertical lines represent lcr-noise of the 
magnetogram signal. 
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Fig. 12. The same as in Fig.HUfor 1600 A flux. 

Table 3. Dependence of the network emission on the magnetic 
flux: exponents of the power-law fits calculated for the images at 
4" and 12" resolution, averaged over 3.5h and 0.5h 



NWflux 


4" resolution 


12" resolution 




3.5h 


0.5h 


3.5h 0.5h 


\B\ - IcallK 


0.6 


0.43 


0.54 0.47 


1^1 ~~ ^1600A 


0.49 


0.53 


0.66 0.44 



(see TableO. Finally, millimeter emission from the network re- 
gions as a function of the absolute value of the magnetic field is 
also best described by a power law (among the considered func- 
tions, see Fig. [T3l . In this case, however, exponents larger than 
unity are obtained and we regard this result with caution and do 
not include it in Table [3] 

The fact that chromospheric emission (be it from Ca n K, the 
1600A range or 3.5 mm continuum) from the internetwork re- 
gions shows almost no dependence on the magnetogram signal 



Fig. 13. The same as in Fig.QTjfor 3.5 mm brightness. 



(see upper panels of Figs. [TT1[T2"1 and [T3l is not really surprising 
given that the magnetogram values in the internetwork do not 
exceed 10 Gauss, with mean of 1.5 Gauss and 1.0 Gauss (for 
0.5h and 3.5h time averages, respectively). Therefore, even after 
creation of 0.5h averaged images, which reduces significantly 
the noise in the magnetograms, only about 3% of IN points dis- 
play B values exceeding the lcr noise level (see left panels of 
Figs. ITTI and [T2l. Even after averaging over 3.5 hours, only 20% 
of IN locations have B values greater than lcr-noise. Data that 
have been so heavily averaged are only of limited value to study 
the internetwork as IN features are commonly characterized by 
much shorter time scales. We also studied the internetwork emis- 
sion using a less strict criterion to discriminate between the cell 
interiors and the cell boundaries. Points with the calcium inten- 
sity not exceeding the value corresponding to the intercept point 
of the two-Gaussian fit to the intensity histogram were selected 
to represent the IN mask. However, the use of this mask does 
not influence the results of the fit. Again, no dependence of IN 
chromospheric emission on the magnetogram signal was found. 
Consequently, the high level of noise in the MDI magnetograms 
used here (compared to the range of the B values in the internet- 
work regions) hinders a reliable study of the relation of internet- 
work chromospheric emission on magnetic field. 

4.6. Chromospheric flux-flux relations 

We studied flux-flux relationships by constructing scatter plots 
for chromospheric intensities. In this analysis excess fluxes for 
K-line and 1600 A emission were considered, employing the 
basal fluxes derived earlier (see Sect. 4.3 and Fig.|7|. For the mm 
intensity we included the zero flux as a parameter in the fitting 
power-law function. Next, using the x 1 criterion we searched for 
the best linear and power law fits to these relationships. The scat- 
ter plots of intensity at 1600 A vs. Can K-line intensity as well 
as Can K-line and 1600 A vs. 3.5 mm intensity are depicted in 
Fig. [14] together with the best power-law and linear fits. For all 
three dependencies the curves depicting both fits are hardly dis- 
tinguishable and the fits are of similar quality (see Fig. FBI. 

For the relation between K-line intensity and intensity at 
1600 A the power-law exponent is 1, which implies that the 
relationship between the excess fluxes is close to linear at a reso- 
lution of 4". The relationship between mm and calcium, as well 
as between mm and UV emissions at the resolution of the mm 
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Fig. 14. Scatter plots of chromospheric fluxes and their power- 
law fits (plotted as solid lines). For comparison we also plot lin- 
ear fits as dotted lines, a) 1600 A excess flux vs. Call K-line 
excess flux, b) Call K-line excess flux vs. 3.5 mm flux, and c) 
1600 A excess flux vs. 3.5 mm flux. 

images, are both best described by power laws with the best-fit 
exponents of 0.7, but as mentioned this exponent is rather un- 
certain since a linear function also represents the data well. 

5. DISCUSSION AND CONCLUSIONS 

We have considered chromospheric disk-center, quiet-Sun emis- 
sion at different wavelengths and its dependence on photospheric 
magnetogram signal, including, for the first time, observations 
at 3.5 mm. Maps in Can K core intensity, the TRACE 1600 A 
channel, BIMA 3.5 mm brightness temperature and (unsigned) 
MDI magnetograms reveal clearly the same pattern, although at 
different spatial resolutions. At the resolution of 12" millime- 
ter brightness resembles quiet-Sun time-averaged emission from 
other chromospheric heights, and demonstrates a dependence on 
the underlying magnetic field similar to that of the other chro- 
mospheric diagnostics. 

The tightest relation between the chromospheric emission 
and photospheric magnetic field is displayed by the Can K 
line. At all wavelengths a power law gives a lower x 1 than 
a linear fit, although for the 1600 A radiation the difference 
is marginal . For Can K t his re sult supports the conclusion 
reached by Schriiver et al. (1989), but based on active region 
data, and by lHarvev & White! d 19991) f or both active and quiet- 
Sun magnetic structures. For Can K ISkumanich et al.l d 1975b 
and lNindos & Zirinl d!998l) found a linear dependence to be ade- 
quate. However, the linear dependencies were obtained after ex- 
cluding weak fields (weaker than 25 G and 20 G, respectively) 
and for a magnetogram signal not exceeding 120 G and 50 G, re- 



spectively. In our study we included both weak fields exceeding 
lcr-noise level of the magnetograms (1.5 or 4 G depending on 
the amount of averaging), and the highest values reached in the 
quiet-Sun data (about 230 G at MDI's spatial resolution). This 
suggests that the non-linearity of the relation found in this work 
is influenced mostly by weak chromospheric emission that re- 
veals a strong dependence on magnetic fields. For instance, if 
we set a B threshold of 20 G we cannot distinguish between the 
quality of the linear and power-law fits for all chromospheric 
emissions analyzed here. Our conclusions agree with those pre- 
sented by Rez aei et"aTI (|2007), who found from the POLIS mea- 
surements at 1" resolution a clear non-linear increase of the cal- 
cium core flux with magnetic flux for values B < 100 G and a 
slow increase for higher values. They also demonstrated that the 
power-law exponent depends strongly on the value of the lower 
threshold. 

The difference between the maximum mag netogram signal 
reache d in our data and that in the data of ISkumanich et al.l 
(fl975l) and lNindos & "Zirinl £l998) may indicate either that their 
quiet-Sun regions are even quieter than the one we analyzed, 
or it may reflect differences in spatial resolution. One more dif- 
ference between our analysis and the others is that we consider 
temporally averaged data rather than snapshot images. By using 
time-averaging we reduce significantly the scatter in the rela- 
tionship and the noise in the magnetograms, which lead to more 
reliable power-law fits. By employing snapshot images we get 
more points in the scatter plots but due to higher magnetogram 
noise level the power-law exponent is even lower (and less re- 
liable) than for the time-averaged images. One should also note 
the difference in the analysis techniques between this work and 
the others applied for building the scatter plots. In this work we 
employed pixel-to-pixel comparison of the images. However, the 
results of the fits are similar if we consider the binned chromo- 
spheric data averaged within a narrow range of magnetogram 
values. On the whole, from the data analyzed here we are able to 
confidently determine non-linear behaviour of 3 chromospheric 
magnetic -intensity relationships. 

The exact power law exponents of the relations between the 
calcium core emission and the photospheric magnetic field dif- 
fers from that found bv lSchriiver et alJd 19891) (0.3 vs. 0.6) . This 
probably has to do with the smaller < B > range available in 
the quiet Sun, although it cannot be ruled out that the Ca n K 
intensity does n ot increase so rapi dly in the quiet Sun as in ac- 
tive regions (see lSolanki et al.|[l991l) . To compare these explana- 
tions, we studied the dependence of the calcium emission on the 
magnetic flux in an active region close to disk center observed 
by the BBSO and MDI on 11.03.2004. The exponent obtained 
from a power-law fit for the AR emission is surprisingly low 
and similar to the QS values, being in the range 0.28-0.32, with 
the exact value depending somewhat on the criterion chosen to 
identify the AR pixels on the map . The use of a pixel s ize and a 
noise threshold similar to that of Schr iiver et al.l dl989l) does not 
increase the value of the power-la w ex ponent and it stil l differ s 
from that of lSchriiver et all dl989l) and (Harvey & White! dl999l) . 
The difference may have to do with the spectral range covered 
by the different observations. 

The reduction in the power-law exp onent as compared to the 
value derived by Schriiver et al. (1989) may be to a large ex- 
tent the result of the inclusion of weak magnetic flux regions, 
which are systematically characterized by magnetogram field 
strengths that are smaller than the true values due to the pres- 
ence of weak neutral lines with null line-of-sight component. 
Unfortunately, we cannot exclude the influence of these hori- 
zontal magnetic fields, which are missed by MDI observations. 
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Such fields, partly se en in the internetwork regi ons by Hinode 
(e.g.lLites etal.l2 007). are very likely small loops (ICenteno et al.l 
2007; Martinez Gonzalez et al. 20 071) and may be due to a loca l 
dynamo dVogler & Schusslerf l2007t Ischiissler & Voglerfl2008l) . 
It is unclear if such low-lying loops really influence the struc- 
ture of the higher-lying chromosphere. However, if they do con- 
tribute to the heating of those layers then this could explain the 
high chromospheric emission from many spatial locations as- 
sociated with small magnetic fluxes in the MDI data (Fig. [7]). 
Unresolved opposite magnetic polarities within the same indi- 
vidual MDI pixels have a similar effect. Both these effects to- 
gether could explain the result that the power law exponent is 
lower when very quiet (low flux) pixels are included. Studies at 
higher spatial resolution and sensitivity are needed to test these 
hypotheses. 

On the whole, chromospheric emission displays a rather dif- 
ferent dependence on the network and internetwork magnetic 
field. The power-law exponents in the relation between the Ca n 
K-line intensity and the net work magnetic flux ar e found to be 
similar to thos e obtained by |Schriiver et al.l (1 19891) for an active 
region and by lHarvev & White! (1 19991) for the quiet Sun net- 
work. The same conclusion holds for the relation between the 
UV emission at 1600 A and magnetic flux in the network re- 
gions of the quiet Sun. For chromospheric fluxes at a resolution 
of 4" it was found that emission from the network entirely de- 
fines the trend for magnetic flux densities stronger than 20 G. 
Our resul ts for the calcium core emission remarkably confirm 
those of iRezaei et al.l (120071) . who obtained for the ensemble of 
quiet-Sun points a low value for the power-law exponent of ^0.2 
(for snapshot images), while for the network emission from re- 
gions with B > 20 G the exponent reaches the value of 0.5 1 (and 
0.45 for B> 10 G). 

A considerable scatter about the mean relation between the 
millimeter flux at 3.5 mm and the magnetic flux density per 12" 
aperture makes a confident study of the mean trend difficult. 
Although the data strongly suggests that the dependence of the 
3.5 mm brightness on magnetic field is non-linear and best de- 
scribed by a power law (among the functional forms considered 
here), but the sizeable scatter does not allow the exponents to be 
determined reliably. The best fit exponent depends strongly on 
the range of < B > considered. 

In the internetwork regions almost no dependence of the Ca n 
K-line, UV at 1600 A and 3.5 mm intensity on the magnetogram 
signal was found. Unfortunately, due to the very low B values 
in the internetwork, only a few percent of IN structures are as- 
sociated with magnetogram signal exceeding the noise level. 
However, it should be noted that analysis of the cell interiors 
and their relation to the underlying magnetic field carried out in 
IRezaei etal |(|2007|)forCa ii H-line with lower noise and higher 
resolution has revealed similar results. This supports the idea 
that the role of magnetic field is negligible for bright internet- 
work grains (e.g. acoustic heating). 

Correlation analysis performed for the quiet-Sun chromo- 
spheric emissions surprisingly revealed a linear relation be- 
tween the C an K-line intensit y and UV brightness at 1600 A. 
According to lRutten et al.l (Il99ll) . who claimed that the relations 
between the fluxes originating in the same atmospheric regimes 
are commonly close to linear, these chromospheric emissions 
are both generated in the chromospheric temperature regime 
from chromospheric ranges of heights. However, solar models 
place UV emission at 1600 A in the temperature minimum re- 
gion, while calcium emission is believed to be generated from 
heights more typical of the lower chromosphere. The power- 



law relations between calcium flux and 3.5 mm flux, as well 
as between 3.5 mm brightness and 1600 A flux, suggest differ- 
ent temperatures of formation for these pairs of emissions, but 
as the power-law exponent do not differ much from unity we 
may assume that these emissions are still close in the heights of 
formation. However, the flux correlation analysis for the consid- 
ered emissions and its conclusions should be treated with cau- 
tion. Due to the very different temperature sensitivity of, on the 
one hand, UV and calcium intensity, and, on the other hand, 3.5 
mm brightness, they represent different contributions from small 
(unresolved) hot magnetic features and cool weak-field regions 
within the same resolution element. In other words, in an inho- 
mogeneous atmosphere the UV samples mainly the hot features, 
while the mm is composed of radiation from both. This could ex- 
plain some of the difference in the behaviour of 1600 A and Ca n 
K excess emission compared to mm emission, and, in particular, 
the large scatter in the flux-flux plots (Fig.[T4"ii. 

Summarizing, clear spatial correlations found between pho- 
tospheric magnetograms, the Can K-line, UV intensity at 
1600 A and 3.5 mm emission indicate that heating in the quiet- 
Sun lower and middle chromosphere maps out the underlying 
photospheric magnetic field rather well. Our results imply that at 
the chromospheric heights covered by the wavelengths investi- 
gated here we deal with the same heating mechanism. However, 
any mechanism that explains chromospheric heating must ex- 
plain the relationship between chromospheric intensity excess 
and magnetic fields. In particular, it should validate both the 
dominating role of magnetic field in the chromospheric emis- 
sion of the network and the absence of correlation between the 
internetwork emission and magnetic flux. 

Of particular interest is the possibility of studying the middle 
chromosphere at the heights where mm emission originates, but 
because of very limited spatial resolution and insufficient sensi- 
tivity a detailed study like that cannot be conducted at present. In 
general, present data on both chromospheric emission and mag- 
netic field are not sufficient to establish physical mechanisms 
that are acting at chromospheric heights leading to the observed 
brightness structures. In this respect, an extraordinarily power- 
ful tool to study the thermal structure of the solar atmosphere 
from the temperature minimum region to the middle chromo- 
sphere will be provided by the Combined Array for Research 
in Millimeter-wave Astronomy (CARMA) and particularly by 
A tacama Large Millimeter Array (ALMA) as was demonstrated 
in lLoukitcheva et al.l d2008l) . 
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